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Abstract The technique of Blind Iterative De-convolution (BID) was used to 
remove the atmospherically induced point spread function (PSF) from short 
Qj ' exposure images of two binary stars, HR 5138 and HR 5747 obtained at the 

^^ ■ cassegrain focus of the 2.34 meter Vainu Bappu Telescope(VBT), situated at 

Vainu Bappu Observatory (VBO), Kavalur. The position angles and separations 
of the binary components were seen to be consistent with results of the auto- 
correlation technique, while the Fourier phases of the reconstructed images were 
consistent with published observations of the binary orbits. 
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<^^ ' 1. Introduction 

On 

0^ ' Atmospheric turbulence degrades the images obtained by ground based astro- 

(— | | nomical telescopes. Schemes like speckle interferometry (Labeyrie, 1970), Knox- 

O |- Thomson algorithm (Knox and Thompson, 1974), shift and add (Lynds et al., 

1976), triple correlation (Lohmann et al., 1983) have been successfully employed 
^ ■ to restore the degraded images. All these schemes depend on the statistical 

^ | treatment of a large number of images. Often, it may not be possible to record 

d • a large number of images within the time interval over which the statistics of 

the atmospheric turbulence remains stationary. In such cases, where only few 
images can be used, there are a number of schemes to restore the image us- 
ing some prior information about the image. The maximum entropy method 
(Jaynes, 1982), CLEAN algorithm (Hogbom, 1974), and BID (Ayers and Dainty, 
1988) are examples of such schemes. In this paper, we employ a version of BID 
developed by P. Nisenson (Nisenson, 1991), on degraded images of two binaries, 
HR 5138 HR 5138 (m„=5.57, Am=0.2), and HR 5747 (m„=3.68, Am=1.5) 
obtained at the 2.34 meter VBT at Kavalur. 

2. Observations 

The 2.34 meter VBT has two accessible foci for backend instrumentation — a 
prime focus (f/3.25 beam) and a cassegranian focus (f/13 beam). The latter 
was used for the observations described in this paper. The cassegranian focus 



>< 



has an image scale of 6.7 arcseconds per mm. This was further magnified to 
~1.21 arcseconds per mm, using a Barlow lens arrangement (Saha et al., 1987, 
Chinnappan et al., 1991). This enlarged image was recorded through a 5 nm 
filter centred on Ha using an EEV intensified CCD camera which provides a 
standard CCIR video output of the recorded scene. The interface between the 
intensifier screen and the CCD chip is a fibre-optic bundle which reduces the 
image size by a factor of 1.69. A DT-2851 frame grabber card digitises the 
video signal. This digitiser resamples the pixels of each row (385 CCD columns 
to 512 digitized samples) and introduces a net reduction in the row direction of 
a factor of 1.27. 

The frame grabber can store upto two images on its onboard memory. These 
images are then written onto the hard disc of a personal computer. The images 
were stored on floppy diskettes and later analysed on a pentium. The observing 
conditions were fair with an average seeing of ~2 arcseconds during the nights 
of 16/17 March 1990. 

3. Data Processing 

The blind iterative deconvolution technique is described in detail in the litera- 
ture (Bates and McDonnell, 1986). Essentially, it consists of using very limited 
information about the image, like positivity and image size, to iteratively arrive 
at a deconvolved image of the object, starting from a blind guess of either the 
object or both the convolving function. The implementation of the particular 
version of BID used by us is described in Nisenson (1991). 

The algorithm has the degraded image c(x, y) as the operand. An initial es- 
timate of the point spread function (PSF) p(x, y) has to be provided. The 
degraded image is deconvolved from the guess PSF by Wiener filtering, which 
is an operation of multiplying a suitable Wiener filter (constructed from the 
Fourier transform P(u, v) of the PSF) with the Fourier transform C(u, v) of the 
degraded image as follows 

0(U,V) = C{U, V) p( Ui „)p.( tti „) + jv-(u, t ,)JV*(u,t)) 

where O is the Fourier transform of the De-convolved image and N is the 
noise spectrum. 

This result O is transformed to image space, the negatives in the image are set to 
zero, and the positives outside a prescribed domain (called object support) are 
set to zero. The average of negative intensities within the support are subtracted 
from all pixels. The process is repeated until the negative intensities decrease 
below the noise. 

A new estimate of the psf is next obtained by Wiener filtering the original 
image c(x, y) with a filter constructed from the constrained object o{x, y). This 



completes one iteration. This entire process is repeated until the derived values 
of o(x, y) and p(x, y) converge to sensible solutions. 

4. Results 

The results for HR 5138 and HR 5747 were arrived at after 350 iterations. Since 
the intensity of the stars were different, the value of the Weiner filter parame- 
ters also had to be chosen accordingly. Figures 1 (a) , 1 (b) , 1 (c) show the speckle 
image, PSF, and deconvolved image of HR 5138 respectively. Figures 2(a), 2(b), 
2(c) are for HR 5747 respectively The companions of HR 5138 and HR 5747 
are separated by 0.27 and 0.20 arcseconds respectively. Measurements of po- 
sition angle (235 and 110 degrees for HR 5138 and HR 5747 respectively) and 
separation were compatible with the values published in the CHARA catalogue 
(McAlister and Hartkopf, 1988). Horch (1994), too found similar results. The 
magnitude difference for the reconstructed objects were 0.04 and 1.65 respec- 
tively for HR 5138 and HR 5747. Although these values compare quite well with 
those published in the Bright Star Catalogue (Hoffleit and Jaschek, 1982) one 
needs to treat many more objects before one can characterise the photometric 
quality of the reconstructions. 

5. Discussion and Conclusions 



The present scheme of BID has the chief problem of convergence. It is indeed 
an art to decide when to stop the iterations. The results are also vulnerable 
to the choice of various parameters like the support radius, the level of high 
frequency suppression during the Wiener filtering, etc. The availability of prior 
knowledge on the object through autocorrelation of the degraded image was 
found to be very useful for specifying the object support radius. In spite of this 
care taken in the choice of the support radius, the psf for each star contains 
residual signatures of the binary sources. Although suggestions have been made 
for improving the convergence (cf. Jefferies and Christou, 1993), these improved 
algorithms require more than a single speckle frame. 

For the present, it is noteworthy that useful reconstructions are possible using 
single speckle frames. Questions regarding their dynamic range and linearity 
can be answered only after examining a wide range of reconstructions. The 
iterative nature of the algorithm does not lend to an explicit estimation of these 
parameters from a limited sample. New developments in camera electronics 
promise the capability to acquire several images within a short time. This 
will also reduce the level of artifacts. The chief achievement in this paper is 
a demonstration of the scientific potential of BID for resolving bright objects 
acquired using simple apparatus. 
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Figure 2: 2(a), 2(b), 2(c) are for HR 5747 respectively. 
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